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In this paper we analyze the gravitational field of a global monopole in the context of 
f{R) gravity. More precisely, we show that the field equations obtained are expressed 
in terms of F(R) = -^^j^- Since we are dealing with a spherically symmetric system, 
we assume that F{R) is a function of the radial coordinate only. Moreover, adopting 
the weak field approximation, we can provide all components of the metric tensor. A 
comparison with the corresponding results obtained in General Relativity and in the 
Brans-Dicke theory is also made. 
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1. Introduction 

It is well known that several types of topological defe cts m ay have been created 
by the vacuum phase transitions in the early uni verse. These include domain 
walls, cosmic strings and monopoles. Global monopoles are heavy objects formed in 
the phase transition of a system composed by a self-coupling triplet of scalar fields 
4>'^ (a = 1,2,3) whose original global 0(3) symmetry is spontaneously broken to 
U{1). The scalar matter field plays the role of an order parameter which outside 
the monopole's core acquires a non- vanishing value. The spacetime associated with 
global monopoles is characterized by a non-trivial topology observed as a deficit 
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solid angle. The properties of the gravitational field produced by this exotic object 
have been investigated by M. Barriola and A. Vilenkin in the context of the general 

o 

relativity. Later, the gravitational field of a global monopoles were analyzed by A. 
Barros and C. Romero in Ref. |4]in the Brans-Dicke theory of gravity in the context 
of weak field approximation. 

On the other hand, the f{R) theories of gravity have been recently suggested as 
a possib le al ternative to explain the late time cosmic speed-up experienced by our 
uni verse. l^l^ Such theories avoid the Ostrogradski's instability that can otherwise 
prove to be problematic for general higher derivatives theories. 

In this work we analyse the gravitational field of the global monopole in the f{R) 
theories of gravity scenario, handling the field equations at the weak field regime. 
In order to simplify the form of the fields equations, we shall express them in terms 
of the function F{R) = ^^j^, in a similar procedure as adopted in Refs.[9]and [TOl 
Since we are dealing with a spherically symmetric static system, we shall express 
F{R) as a function of the radial coordinates, r, by defining F{R{r)) = J-{r). 



2. The Global Monopole Model in General Relativity 

The simplest model which gives rise to a global monopole is described by the La- 
grangian density below:- 

c = ^d.cj^'^d^^r - iK^^r - ■ (1) 

Coupling this matter field with the Einstein equations, spherically symmetric solu- 
tions for the fields equations can be obtained by adopting for the line element and 
the matter fields the Ansdtze below: 

ds^ = B{r)dt^ - A{r)dr^ - r^{de'^ + sin^ Odip^) (2) 

and 

r^vHr)- (3) 
r 

with x°'x°' = . 

For points far from the monopole's core, the solutions for the functions B{r) 
and A{r) are given by 

B = A-^7iil- SirGrf - 2GM/r , (4) 

being the mass parameter M « Mcore- What concerns the matter field, we have 
h{r) « 1. Moreover, in Ref. [3] the authors have shown that for points outside 
the monopole's core, the energy-momentum tensor associated with a static global 
monopole can be approximated as 

TA'«diag(^,^,0,o) . (5) 
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3. The /(-R) Gravity in the Metric Formalism 

The action associated with the modified tlicories of gravity coupicd with matter 
fieids is given by: 

S^^ ( d^X^gf{R)+Sra , (6) 



2k _ 

where /(i?) is an analytical function of the the Ricci scalar, R, k = SttG and Sm 
corresponds to the action associated with the matter fields. By using the metric 
formalism, the field equations become: 

in which T/;^ is the geometric energy-momentum tensor, namely 



T'^. = [Ig,,^ ifiR) - m)R) + V"V^F(i?) - <?^.g„^)| 



(8) 



The standard minimally coupled energy- momentum tensor, TJ^l, derived from 
the matter action, is related to by 

= TirjF{R) . (9) 
Thus the field equations can be written as 

F{R)R^.u - lfiR)9^.. - V^V,^^(i?) + g^.nF{R) = kTJZ ■ (10) 

Taking the trace of the above equation we get 

F{R)R~2f{R)+3aF{R)^ kT"" , (11) 

which expresses a further scalar degree of freedom that arises in the modified theory. 
Through this equation it is possible to express f{R) in terms of its derivatives and 
the trace of the matter energy-momentum tensor as follows: 

f{R) = I iF{R)R + iOFiR) - kT"') . (12) 
Substituting the above expression into (llOp we obtain 



F{R)R^, - V^V,F(i?) - KfjJl = ^ [F{R)R - DF{R) - «:T'"] . (13) 



From this expression we can see that the combination below 

F(R)R,,^ - V. V„i^(i?) - kT,™ 

with fixed indices, is independent of the corresponding index. So, the following 
relation 

- = , (15) 

holds for all fi and v. 
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4. The Global Monopole in f{R) Theories 

In this section wc shall derive the field equations associated with a global monopole 
system in modified theories of gravity. 

By adopting a spherical symmetry Anstaz for the metric tensor given in ([5]), the 
non- vanishing components of the Ricci tensor read: 



Rq 



Rl 



1 
2 

1 

'2 



/\2 



B'A' 



B" 2B' 



2B^A 2BA^ BA BAr 



B'A' 
2BA^ 



d2 A' 
The spherical symmetry requires that 



2A' 
B' 



B" 
BA 

1 



/^2 



{B') 
2B^A 



2rBA r^A 



d2 r>3 

^2 — -"-3, 



which implies that the scalar curvature will be given by 
2 



R 



^ , B' A'\ IB' „ fB' A' 
r^A\'-''+[B-A)'-4B'[B^A 



1 

2 B 



(16) 
(17) 
(18) 

(19) 

^ \ ■ (20) 



In all the above equations the primes corresponds derivative of the function with 
respect to the radial coordinate. By using the identity ([T5|) along with the energy- 
momentum tensor given by ([5]), we can construct two linearly independent differ- 
ential equations: 







fB' 






'a) 












'a) 



and 



2rT -rT' 



AB + AAB - ArB— + 2r^B'— - r^B' ( — 







(21) 



B' A' 
B + A 



2r^B 



: B' A'\ AABki? 



(22) 



where we have expressed F{R) as F{R{r)) = ^{r). 
Defining 



p. 



B' A' 



we can write the field equations as follows 



^^^_1^ 
r J- 2 J- ^ ' 



(23) 



(24) 
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and 

F' T' 
- 4B + - 4rB— + 2r^B' — 
J" T 

+ 2r^B" - r^B'l3 + 2Br(i - ^^^I^ = o . (25) 

We can verify that if we plug oj — into the equation (9) of Ref. l4j and replace 
4>{r) by J-{r), the equation (f24| is recovered. This emphasizes the equivalence of 
the metric f{R) gravity with a Brans-Dicke gravity possessing a parameter oj — 0, 
as pointed out in Ref. 1141 Furthermore this fact shows the relevance of working with 
the fields equations in terms of F{R) when one wants to compare with the results 
obtained via Brans-Dicke gravity, since the role of the scalar degree of freedom in 
f{R) gravity is played by the function F{R). 

5. Solutions in the Weak Field Regime 

Now let us consider the weak field approximation in the field equations by assuming 
that B{r) = 1 + b{r) and A{r) = 1 + a{r) with \b{r)\ and \a{r)\ smaller than unity. 
Furthermore, let us consider that the modification theory of gravity corresponds 
to a small correction on General Relativity (GR), so that J-{r) = 1 + ip{r), with 
mr)\ « 1 

Adopting these approximations it is possible to verify that up to the first order 
we may write, 

J" l + i/)~'^'J" 1 + V'~^ 
B' b' A' a' 

— = r ~b' , — = « a' . 26 

B l + b ' A 1 + a ^ ' 

Thus the approximated equations for (PT|) and (|22p read 

^ = ^" (27) 
r 

and 

4a - 4rV'' + 2r{a + b') - 2r%" - 4(1 + a + & - ^)k7?^ . (28) 

Specific results for the fields equations can be obtained by adopting for i!{r) a 
specific Ansatz. Considering the simplest analytical function of the radial coordinate, 
namely ^{r) = ^^r, we shall have the following result for Eq. (P7|) : 

P , . 

- = , 29 
r 

which means 



a' + = 
a + = cq , 



(30) 
(31) 
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where cq is an integration constant that we choose as cq = 0. Therefore 

a{r) = -b{r) . (32) 

So (|28p will be written as 

hrH" - 5 - rV'o - (1 - MKif = , (33) 

whose solution is 

h(r) = — + C2r^ — nrf — ilior[l — nif') , (34) 
r 

where ci and C2 are integration constants. Since we are not dealing with a theory 
with a cosmological constant we may set C2 — 0. Furthermore, we must choose 
Ci = —2GM in order to recover the Newtonian potential present in the Barriola 
and Vilenkin solution. Therefore, 

9C AT 

B{r) = 1 - =^ - 8^G?72 - 1^0^(1 - SttGj]^) , (35) 
r 

since k = SnG. For a typical Grand Unified Theory the parameter 77 is of the order 
10^^ Gev. So, STTGif w 10"^.^ This allows us to neglect the term tp^r x SnGi]^ in 
the above result since \ipof \ << 1- Then we may write 

B{r) = 1 - - SttGt]^ - Vor . (36) 



From equation (j29p we have 



Air) = , (37) 

B{r) 



where ao is an integration constant. Rescaling the time coordinate we can set ao = 1. 
Then, 



A{r) 



2GM 2 ■ 

1 sttGi] — tpor 



(38) 



As previously assumed by Barriola and Vilenkin, here we also drop out the mass 
term in (l36t and (|38l) as it is negligibly small at astrophysical scale. Thus, we have 

B{r) = 1 - SnGri^ - V'o^ (39) 

and 

Air) ^ {1 - SnGri^ - Tpory^ ■ (40) 

Using the weak field approximation in (|20l) we can obtain the Ricci scalar curvature 
associated with the solution above: 

= (41) 
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From J-{r) and R{r) one can determine, in principle, F{R) and finally f{R). The 
explicit form obtained for /(i?) is given by 

R 



fiR) = i? - 3V^^ In ( — ) - 2V^o V9^o " 16^Gry2i? 



^/QV'g - m7TGT]^Ra - 3^0 



■ 3V'o In 



(42) 



where i?o is a constant. It is straightforward to verify that if ipo is set to be posi- 
tive, the function above satisfies th e fo llowing stability conditions required for any 
physically relevant /(i?) theoryP^^l^l 

• ^-dW^ > tachyons); 

• > (no ghosts); 

• limn^oo ^ — ^ and Mmji^ao ^ = (GR is recovered at early times), 

where A = A(i?) is defined as A = f{R) - R. 

Now we shall show that the line element described by the obtained functions 
A{r) and B{r) are conformally related to the ordinary global monopole solutions^ 
So let us consider the following coordinate transformation: 

B{r) = p{r*) (1 - SirCrf) , (43) 
A{r)dr^ ^ p{r*) (l + SttG??^) {dr* f , (44) 
r =pi/2(r*)r* , (45) 

where h{r*) is an arbitrary function of r* to be determined, and h{r*) = 1 + q{r*) 
with |<7(r*)| << 1. Differentiating equation ([45|) we have 

dr"" =(l + r*^ + <^dr*^ . (46) 
\ dr* J 

Substituting the above equation into (|44p and keeping only linear terms in q{r*), 
ipor and Gr]^, we obtain the following result for q{r*): 

q{r*) = -ijjor* , (47) 

then 

p{r*) ^ 1 - t/jor* . (48) 
Thus wc can write the line element ([2]) in the coordinate r* as follows: 
ds^ = (1 - t/jor*) [(1 - SttGt]^) dt^ - (l + SirGrj"^) dr*^ 



r*2 (de^+sin'^ 



edcj)^)] . (49) 



^An analogue procedure has been adopted in Ref. |4l 
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Rescaling the time coordinate and redefining the radial coordinate as r = 
(l + AttGi]'^) r* , we arrive at the line element below: 

ds^ tpor) [dt^ - dr^ - (l - SttGt/^) (d0^ + sin^ ed(P^)] . (50) 

An important feature arises if we are interested in analyzing the deflection of 
light in this metric. As it is well known, the deflection angles are always preserved 
for two metrics related by a conformal transformation. Therefore, the deflection of 
a light ray by the monopole in the present modified gravity will be the same of that 
one previously obtained in Refs. |3l 

5<l)^8TTGr]^l{d + l)-\ (51) 

where d and I are the distances from the monopole to the observer and to the source, 
respectively. In this analysis it was assumed that the light ray trajectory lies at the 
equatorial plane 9 = ^. 

6. Conclusions 

In this work we have analyzed the gravitational field of a global monopole in the 
f{R) gravity scenario in the metric formalism. In this formalism, this modified 
gravity contains a massive scalar degree of freedom in addition to the familiar 
massless graviton and it turns out to be equivalent to a Brans-Dicke theory. In order 
to simplify our analysis we have considered solutions in the weak field approximation 
which implied that we are considering a theory which is as a small correction on 
GR. The latter condition was explicitly considered by assuming T{r) = 1 + ipor, for 
\^or\ « 1. 

Following the above mentioned approach, the solutions found by us correspond to 
small corrections on 500 and gn components of the metric tensor, only. Being these 
new components given in p6|) and (|38|) , respectively. From the results obtained, we 
can observe that also, a small correction on the Ricci scalar, given now by (|4ip. 
takes place. Moreover, we have also verified that these solutions are conformally 
related to that one previously obtained by M. Barriola and A. Vilenkin in Ref. [3|as 
shown in ([50| , what ensures us that the defiection of light in these two spacetimes 
will be the same. 

Before finishing this section, we would like to point out that by using the weak 
field approximation, we were able to provide an explicit expression for the function 
f{R), given in (011). The non- linear function f{R) is not a problem if we consider 
low curvatures (i.e., late time in the matter dominated era). Also, by assuming that 
ipo is a positive parameter, we can affirm that f{R) fulfills two important stability 
conditions, namely, non-tachyons and non-ghosts requirements. 
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